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1.  MAGNETIC  STORM  VARIATIONS  IN  THERMOSPHERIC  DENSITIES, 

TEMPERATURES,  AND  WINDS 

At  the  NOAA  "Workshop  on  Satellite  Drag,"  March  18-19, 
1982,  Boulder,  Colorado,  a  paper  entitled  "Geomagnetic  Storm 
Variations  and  Prediction  of  Low-Perigee  Satellite 
Ephemerides"  by  J.M.  Forbes  was  presented.  A  written  version 
has  appeared  in  the  proceedings  of  the  workshop.  In  this 
paper,  the  spatial/temporal  characteristics  of  density 
perturbations  in  the  lower  thermosphere  (150-250Km)  during 
geomagnetic  disturbances  are  discussed  within  the  context  of 
the  low-perigee  satellite  ephemeris  prediction  problem.  The 
variability  of  atmospheric  density  (satellite  drag)  is 
significantly  enhanced  at  high  latitudes  due  to  the  proximity 
of  magnetosphere/ionosphere  source  mechanisms  for  the  neutral 
atmosphere  disturbances.  In  addition,  for  precise  orbital 
computations  winds  in  excess  of  about  300msec_1(which 
commonly  occur  at  high  latitudes)  must  also  be  considered  as 
part  of  the  drag  effect.  It  is  demostrated  that  precise 
low-perigee  ephemeris  predictions  require  accurate 
forecasting  of  the  phase  as  well  as  spatial  structures  of 
atmospheric  density  and  wind  disturbances.  It  was 
recommended  that  future  advances  in  predictive  capability  of 
atmospheric  drag  necessitate  (1)  Comprehensive  numerical 
modelling  efforts  which  tie  together  all  of  the  relevant 
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high-latitude  processes  in  a  self-consistent  manner,  and 

which  can  be  utilized  to  provide  insight  into  the  development 

of  an  improved  magnetic  index  of  density  variations,  and 

guide  the  development  of  semi-empirical  models  for  predicting  » 

stellite  drag  variations;  and  (2)  Concerted  experimental 

efforts  to  obtain  simultaneous  drag  and  electrodynamic  data 

(incoherent  scatter  measurements  of  temperature  and  drifts; 

satellite  measurements  of  particles,  fields  and  currents;  and 

ground  magnetometer  measurements). 

Dr.  Forbes  and  Mr.  J.  Slowey  attended  a  Technical 
Interface  Meeting  at  AFSCF ,  Sunnyvale,  CA  on  26-27  April 

1982,  which  was  held  for  the  purpose  of  coordinating  AFSCF' s 
Data  System  Modernization  <DSM)  efforts  with  its  prime 

contractor,  IBM.  A  point  of  discussion  concerned  the 

Jacchia-Roberts  analytic  model,  which  was  originally  fitted 
to  the  J70  model,  whereas  a  second  version  of  the  model  was 
developed  by  IBM  which  was  fitted  to  the  J71  model.  J.  Slowey 
recommended  that  the  version  of  the  Jacchia-Roberts  model 
fitted  to  the  J70  model  be  retained  in  the  system,  but  that 
equations  representing  low  altitude  geomagnetic  activity  in 
the  J71  be  substituted  for  those  in  the  J70  model.  Dr. 

Forbes  gave  some  reasons  why  it  would  not  be  advisable  to  add 
the  MSIS-77  model  to  their  system  at  this  time.  These 

included:  * 


VV.  .V.  -V-V  V  V  W-- 


(a)  The  MSIS  represents  a  fit  to  composition  and  temperature 
measurements,  not  total  mass  densities; 

(b)  Data  utilized  to  construct  MSIS  covers  a  limited  range 
of  solar  activity  and  altitude,  as  compared  to  the 
Jacchia  models; 

(c)  The  geomagnetic  effect  in  MSIS-77  is  parameterized  by 
the  daily  Ap  index,  which  is  incapable  of  reflecting 
time  variability  as  well  as  the  3-hour  Kp  index  utilized 
in  the  Jacchia  models; 

(d)  MSIS-77  requires  about  25%  more  CPU  time  that  J71;  and 

(e)  The  Millstone  Hill  (MH)  temperatures  used  in  MSIS-77 
originated  from  the  1-pulse  experiment  at  MH,  which 
contain  a  50K  bias  (based  on  more  accurate  2-pulse  data 
analyzed  since  development  of  MSIS-77.). 

Dr.  Forbes  also  presented  comparisons  of  prediction  errors 
for  various  models  and  satellite  orbits  which  have  been 
examined  over  the  years  within  AFGL.  In  particular,  the 
Harris-Priester  Model,  which  was  discussed  for  possible 
consideration  by  SCF  and  IBM  as  part  of  DSM  was  shown  to 
exhibit  considerably  larger  prediction  errors  and  about  twice 
the  CPU  time  than  J71. 

Dr.  Forbes  and  A.  Kantor  of  AFGL  attended  a  meeting  of 
the  DARPA  Committee  for  Autonomous  Positioning  of  GPS  at  Hq 
ONR  on  May  1982.  At  this  meeting  possible  satellites,  in 


» 


addition  to  GPS,  for  which  the  SHAD  (Stellar  Horizon 
Atmospheric  Dispersion)  navigation  system  could  be  used  were 
discussed.  Some  of  these  satellites  are  classified. 
Density  variability  in  the  20-50Km  altitude  region  could 
provide  a  limitation  to  the  navigation  accuracy.  A.  Kantor 
discussed  his  planned  study  of  day-to-day  and  seasonal 
variabilities,  including  the  statistical  quantities  to  be 
provided,  and  Dr.  Forbes  presented  his  plans  to  provide 
information  on  the  tidal  variations.  In  addition,  in 
response  to  Dr.  F.  Quelle' s  (ONR)  request  for  possible 
science  that  could  be  derived  from  the  SHAD  experiment,  Dr. 
Forbes  offered  a  number  of  suggestions,  falling  within  the 
following  scientific  areas: 

(a)  Global  Climatology 

(b)  Polar  Stratospheric  Warmings 

(c)  Planetary-Scale  Waves 

(d)  Equatorial  Waves 

(e)  Mesospheric  Particulate  Layers 

(f)  Minor  Constituents 

which  were  subsequently  incorporated  with  other 
recommendations  into  a  letter  from  Dr.  Champion  to  Dr. 
Quelle . 

During  1982  the  capability  was  developed  to  derive 
neutral  thermospheric  temperatures  from  ionospheric 
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parameters  (Te,  Ti,  Ne)  measured  by  incoherent  scatter  radars 
at  Millstone  Hill  .(42°)  and  Arecibo  (18°H).  The  purpose 
of  this  project  was  to  examine  the  equatorward  propagation 
and  extension  of  high-latitude  geomagnetic  disturbance 
effects.  To  develop  the  above  capability  involved  the 
following  tasks: 

(1)  Acquisition  of  the  following  magnetic  tapes  of  Arecibo 
data  from  the  World  Data  Center: 

(a)  19/129  (August  1974  -  May  1977) 

(b)  N35  (December  1971  -  December  1972) 

(c)  N12  (July  1966  -  December  1970) 

The  tapes  initially  acquired  contained  formatting  errors 
and  additional  problems.  After  extensive  attempts  at 
analysis  and  conversations  with  the  World  Data  Center, 
reformatted  tapes  were  acquired  and  subsequently 
analyzed. 

(2)  Development  of  programs  to  read  and  plot  Te,  Ti,  and  Ne 

off  acquired  tapes.  A  sample  plot  is  illustrated  in 
Figure  1.  For  every  set  of  vertical  profile 

measurements  (available  on  the  order  of  every  15 
minutes)  of  Te,  Ti  and  Ne,  a  determination  had  to  be 
made  of  the  quality  and  acceptability  of  the  data  for 
further  analysis.  Sometimes  instrument  or  power 
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TIME  =  64.80 


TIME  *  65. 17 


Figure  1.  Sample  vertical  profile*  of  ionospheric  parameters 
derived  from  Arecibo  incoherent  scatter 
measurements  on  day  260  of  1974. 


problems  suggested  the  rejection  of  data.  Occasionally, 
unrealistic  or  inexplicable  observations  (for  instance  Ti  > 
Te)  had  to  be  rejected.  The  editing  of  Te,  Ti  and  Ne  data 
prior  to  neutral  temperature  determinations  (see  following 
paragraph)  is  a  tedious  step  in  the  analysis. 

(3)  Acquisition  of  NEUTEMP  program  from  Dr.  W.  Oliver  of 
Millstone  Hill  and  program  redevelopment  for  use  on  the 
AFGL  CDC  machine  for  calculation  of  Tn.  The  Millstone 
Hill  facility  uses  a  Harris  computer.  The  NEUTEMP 
program  iteratively  solves  the  ion  thermal  balance 
equation  (given  Te,  Ti,  Ne)  to  produce  neutral 
thermospheric  (exospheric)  temperatures. 

(4)  Tailoring  of  the  Multiple  Linear  Regression  Analysis 
(MLRA)  Program  and  development  of  plotting  package  for 
Millstone  Hill  neutral  exospheric  temperature  analysis 
on  the  CDC  machine. 

The  first  Arecibo  experiment  analyzed,  and  the  one  used  to 
develop  and  test  the  above  programs,  is  day  260  of  1974. 
This  day  was  characterized  by  an  F10.7  value  of  89  and  a 
daily  magnetic  Ap  index  of  3  (Kp  less  that  2-  throughout  the 
day).  The  neutral  exospheric  temperatures  obtained  are 
illustrated  in  the  lefthand  portion  of  Figure  2.  The  MLRA 
yielded  the  following  Fourier  components:  mean  =  842K, 
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diurnal  *  158K,  semidiurnal  *  24K,  and  terdiurnal  *  22K. 
Local  times  of  maxima  were:  diurnal  *  (15.2),  semidiurnal  * 
(3.0,  15.0)  and  terdiurnal  *  (4.0,  12.0,  20.0).  The 
following  day  (September  18)  was  characterized  by 
significant  geomagnetic  activity,  with  Kp  values  of  up  to 
5+ .  The  quiet-day  least-squares  fit  was  extrapolated  into 
the  second  day  to  delineate  the  different  temperature 
response  during  the  magnetically  disturbed  period  (See  Figure 
2).  Temperature  differences  from  the  quiet-day  behavior  in 
the  range  of  100-200K  are  clearly  evident. 

Exospheric  temperatures  were  also  derived  from 
incoherent  scatter  data  at  Millstone  Hill  (42°N)  during  the 
geomagnet ically  disturbed  (Kp  <  5o)  day  of  June  3,  1971.  A 
magnetically  quiet  day  (Kp  <  2-;  F10.7  *  108.3)  with  similar 
F10.7  was  also  analyzed  as  a  reference  day.  Unfortunately, 
the  quiet-day  data  contained  strong  oscillations  suggesting 
the  presence  of  gravity  waves,  rendering  questionable  the  use 
of  a  fit  to  these  data  as  a  quiet-day  reference.  Therefore, 
the  separation  of  the  magnetic  storm  contribution  from  the 
June  3,  1971  data  was  not  possible. 

The  next  step  in  the  incoherent  scatter  project  was  to 
have  been  the  analysis  of  several  storms  simultaneously 
observed  at  Millstone  Hill  and  Arecibo  during  the  solar 
maximum  period  1979-1981.  However,  besides  the  March  22, 
1979,  storm  to  be  described  below  in  conjunction  with  the 


analysis  of  satellite  accelerometer-measured  densities, 
this  effort  was  never  completed.  The  reasons  for  this  were 


twofold.  First,  initial  postponement  of  the  effort  occurred 
due  to  the  higher  priority  given  to  the  interpretation  of 
satellite  winds  and  densities  during  geomagnetically  active 
periods,  as  described  below.  Then,  approximately 
simultaneous  with  the  principal  investigator  moving  to  Boston 
University,  a  new  and  unanticipated  effort  relating  to  the 
analysis  of  data  and  development  of  inputs  to  the  Air  Force 
Reference  Atmosphere  in  the  80-120Km  height  region  was  given 
first  priority  by  the  Air  Force.  The  subsequent  subcontract 
to  Dr.  Forbes  at  Boston  University  only  provided  for  (a)  the 
80-120Km  region  data  analysis  and  (b)  satellite  density  and 
wind  analysis  during  magnetically  disturbed  conditions,  and 
did  not  fund  completion  of  the  incoherent  scatter  work.  Work 
completed  under  item  (b)  is  reported  below  and  the  results  of 
(a)  are  included  in  Section  3. 

The  work  involving  analysis  and  interpretation  of 
satellite  accelerometer-measured  densities  and  winds  under 
geomagnetically-  disturbed  conditions  is  reported  fully  in 
the  following  publications: 

Marcos,  F.A.,  and  J.M.  Forbes,  Thermospheric  Winds  from 
the  Satellite  Electrostatic  Triaxial  Accelerometer  System,  J. 
Geophys.  Res.,  90,  6543-6552,  1985. 


in 


Forbes,  J.M.,  and  F. A.  Marcos,  Thermospheric  winds, 
Densities,  and  Temperatures  during  an  Isolated  Magnetic 
Storm.  (CDAW-6  Interval),  J.  Geophys.  Res.,  accepted  for 

publication,  manuscript  in  revision,  1985. 

In  the  first  paper  a  new  thermospheric  wind  measurement 
technique  is  reported  which  is  based  on  a  Satellite 

Electrostatic  Triaxial  Accelerometer  (SETA)  System  capable 
of  accurately  measuring  accelerations  in  the  satellite's 
in-track,  cross-track,  and  radial  directions.  Winds 
measured  between  170  and  210  Km  during  a  5-day  period  of 
mostly  high  geomagnetic  activity  are  analyzed  to  demonstrate 
the  potential  contributions  of  SETA  data  to  studies  of 
thermospheric  dynamics.  The  data  are  consistent  with  a 
two-cell  polar  circulation  pattern  characterized  by  a 
trans-polar  flow  parallel  to  the  1600  h/0400  h  geomagnetic 
local  time  meridian,  and  return  flows  in  the  late  morning  and 
late  evening  sectors.  The  flow  pattern  was  found  to  be 

asymmetric  in  that  it  is  displaced  about  5-10°  latitude 
towards  the  noon  (MLT )  sector,  and  the  evening  cell  is 
somewhat  more  diffuse  than  the  morning  cell.  The  system 
also  covers  a  greater  area  of  the  polar  cap  and  is  more 

intense  during  active  (Kp  >  5°)  than  quiet  (Kp  3°) 
geomagnetic  conditions.  Average  trans-polar  flow  velocities 
are  characteristically  of  order  150  ♦  75  ms_1for  Kp  —  2° 
and  375  ♦  100  ms-lfor  Kp  _  5Of  the  stated  variabilities 


representing  1<T  deviations.  Details  may  be  found  in  the 
above-cited  publication. 

In  the  second  work,  thermospheric  winds  and  densities 
form  SETA  (Satellite  Electrostatic  Triaxial  Accelerometer) 
measurements  and  exospheric  temperatures  form  Millstone  Hill 
incoherent  scatter  radar  data  are  examined  during  an  isolated 
magnetic  storm  occurring  on  March  22,  1979,  a  period  which 
coincides  with  the  CDAW-6  interval.  A  polar  thermospheric 
wind  circulation  suggestive  of  a  two-cell  horizontal 
convection  pattern  is  reflected  in  the  cross-track 
acceleration  measurements.  Winds  are  generally  trans-polar 
and  appear  to  flow  parallel  to  the  1600h/0400  MLT  meridian 
with  return  flows  in  the  morning  and  evening  sectors  near 
57.5°  +  7.5°  and  70°  +  7.5°  geomagnetic  latitude, 
respectively.  Winds  are  of  order  100-200  ms_1control 
during  quiet  periods,  and  attain  maximum  speeds  between  300 
and  600  ms_1during  the  storm.  The  two-cell  pattern  is 
distinctively  more  well-ordered  in  geomagnetic  rather  than 
geographic  coordinates.  The  time  delay  between  maximum 
magnetic  disturbance  as  reflected  by  the  3-hour  ap  index  and 
full  set-up  of  the  circulation  pattern  is  somewhere  between 
0-3  hours,  the  uncertainty  being  due  to  the  discrete  nature 
of  ap  and  the  90-minute  orbital  period.  The  circulation 
pattern  persists  almost  unattenuated  for  about  6  hours  after 


the  magnetic  disturbance  has  returned  to  quiet  levels, 
reflecting  the  so-called  'flywheel'  effect. 

The  density  response  is  highly  asymmetric  with  respect 
to  its  day/night  behavior.  At  high  geomagnetic  latitudes  (> 
60°)  the  daytime  density  variation,  about  a  40%  increase 
form  quiet  levels,  distinctively  reflects  the  response  to  an 
increase  in  magnetic  activity  with  a  time  delay  of  less  than 
3  hours.  At  lover  latitudes  the  response  is  smaller  — 20%) 
and  less  well-defined,  occuring  with  time  delay  of  about  6  + 
2  hours  near  10-20°  latitude.  The  nighttime  density 
response,  on  the  other  hand,  is  not  so  veil  defined  poleward 
as  it  is  equatorvard  of  60°  latitude.  The  time  delay 
increases  from  about  2+2  hours  to  6  +  2  hours  from  high  to 
low  latitudes. 

An  exospheric  temperature  response  occurs  at  Millstone 
Hill  (42°N)  with  an  amplitude  of  210  K  and  time  delay  of 
between  0  and  2.5  hours. 

Latitude  structures  of  the  density  response  at 
successive  times  following  the  substorm  peak  suggest  the 
equatorvard  propagation  of  a  disturbance  with  a  phase  speed 
of  between  300  and  600  ms-1.  A  deep  depression  in  the 
density  at  high  latitudes  (>70°)  is  evident  in  conjunction 
with  this  phenomenon.  The  more  efficient  propagation  of 
the  disturbance  to  lower  latitudes  during  the  night  may  be 
connected  with  the  'midnight  surge'  effect,  a  reinforcement 
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of  southward  winds  associated  with  high  latitude  heating  and 


solar  EUV  heating,  as  opposed  to  cancellations  between  these 


meridional  flows  during  the  daytime.  Details  of  this  study 


may  be  found  in  the  above-cited  manuscript. 


2.  ATMOSPHERIC  TIDES  STUDIES 


Atmospheric  tides  constitute  much  of  the  meteorology  and 
variability  in  atmospheric  fields  above  70  Km  altitude,  and 


as  such  represent  a  major  uncertainty  in  specifying 


environmental  conditions  for  a  variety  of  Air  Force 


missions.  Experimental  data  often  suffer  from  inadequate 


coverage  in  local  time  and/or  latitude  to  properly  delineate 


tidal  motions,  so  greater  reliance  on  theory  and  numerical 


simulations  must  be  made  in  this  discipline  of  study.  The 


following  work  was  performed  to  better  understand  the  nature 


of  tidal  oscillations  and  the  uncertainty  which  they 


contribute  to  density  and  temperature  specifications  in  the 


upper  atmosphere. 


A  paper  entitled  "Neutral  Temperatures  from  Thompson 


Scatter  Measurements:  Comparisons  with  the  CIRA  (1972)"  by 


J.  M.  Forbes,  M.  E.  Hagan,  and  K.  S.  W.  Champion  was  prepared 


for  presentation  at  the  COSPAR  Meeting,  Ottowa,  Canada,  May, 


1982,  and  was  published  in  Space  Research,  Volume  3,  1983. 


In  this  paper.  Neutral  exospheric  and  lower  thermospheric 


(100-130  km)  temperatures  from  Thomson  scatter  measurements 
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at  Millstone  Hill  (42°N)  are  compared  with  CIRA 
temperatures  with  a  view  towards  identifying  deficiencies  in 
the  CIRA  and  recommending  revisions.  CIRA  is  found  to  model 
the  observed  diurnal  mean  temperatures  (T0)  to  within  10% 
over  a  wide  range  of  solar  conditions  (75<  F10.7<  250),  but 
consistently  underestimates  the  diurnal  temperatures  with 
maximum  deviations  approaching  50%  of  observed  amplitudes 
(180-240  K)  at  solar  maximum  (1200  K  <  T0  <  1400  K) .  The 
observed  semidiurnal  amplitudes,  which  lie  in  the  range  of 
20K-80K,  are  always  underestimated  and  frequently  by  more 
than  50%.  In  the  lower  thermosphere,  tidal  oscillations  of 
temperature  of  order  20K-40K  occur  which  are  not  modelled  by 
CIRA.  In  addition,  an  analysis  of  exospheric  temperatures  at 
Millstone  Hill  during  a  magnetic  disturbance  indicates  a 
response  within  1-2  hours  from  storm  onset,  whereas  CIRA 
assumes  a  6.7  hour  delay.  As  a  result,  and  although  some  of 
these  deficiences  are  addressed  by  the  more  recent  MSIS 
model,  there  exists  a  sufficient  data  base  to  recommend 
several  additional  revisions  to  the  CIRA  temperatures  at  this 
time. 

A  review  of  thermosphere  tides  was  prepared  as  AFGL 
Report  TR-82-0264  entitled  "Tides  in  the  Thermosphere:  A 
Review"  by  J.M.  Forbes  and  K.S.W.  Champion.  In  it  a 
comprehensive  review  of  recent  theoretical  and  observational 
accomplishments  relating  to  diurnal  and  semidiurnal  tidal 


oscillations  of  neutral  winds,  temperature,  density,  and 
composition  above  100  Km  is  presented.  Topics  emphasized 
include:  (1)  Recent  theoretical  studies;  (2)  Solar  cycle, 
seasonal,  and  latitudinal  variations  in  tidal  oscillations  of 
temperature  and  winds  as  inferred  from  Thomson  scatter 
measurements;  (3)  Tidal  variations  in  total  mass  density  and 
composition  as  inferred  from  satellite  accelerometer  and  mass 
spectrometer  measurements;  (4)  Comparison  of  recent 
theoretical  models  with  the  above  observations;  (5)  The 
relative  influence  of  in-situ  and  propagating  tides  in 
determining  the  total  semidiurnal  thermospheric  tide;  and  (6) 
Propagating  tides  of  lower  atmosphere  origin  as  a  source  of 
mean  momentum  and  heat  in  the  lower  thermosphere. 

As  an  input  to  the  SHAD  (Stellar  Horizon  Atmospheric 
Dispersion)  Program,  a  draft  report  was  written  by  Dr. 
Forbes  entitled  "Density  Variability  due  to  Atmospheric  Tides 
below  140  Km".  In  this  report  equations  and  a  computer 
algorithm  were  derived  for  the  tidal  oscillations  in  density 
compatible  with  temperatures  from  the  numerical  tidal  model 
of  Forbes.  Vertical  profiles  and  tables  of  tidal  amplitude 
and  phase  corresponding  to  0°,  30°,  and  60°  latitude 
under  equinox  conditions  were  included  for  the  diurnal  and 
semidiurnal  components,  respectively.  Diurnal  amplitudes 
were  found  to  be  of  order  1-2%  in  the  stratopause  region 
(ca.  50  km)  and  5-15%  in  the  mesopause  region,  attaining 
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amplitudes  of  7-25%  near  110  Km  where  the  first  diurnal 
propagating  (1,1)  mode  begins  to  experience  severe  molecular 
dissipation.  The  phase  structure  at  0°  and  30°  latitude 
in  fact  revealed  a  nominal  25-30  Km  vertical  wavelength 
characteristic  of  the  (1,1)  mode.  A  constancy  of  phase  with 
height  above  50  Km  for  latitudes  >~  30°  was  noted  to  be 
indicative  of  the  trapped  (1,-2)  tidal  mode. 

The  semidiurnal  oscillations  were  found  to  be 
significantly  smaller  than  the  diurnal  component  below  100 
Km,  of  order  .2-1%  around  the  stratopause  and  2-6%  in  the 
mesopause  region.  Amplitudes  of  20-25%  are  attained  in  the 
vicinity  of  120  Km.  Semidiurnal  vertical  wavelengths  below 
70  Km  are  much  larger  than  for  the  diurnal  tide, 
characteristic  of  the  fundamental  (2,2)  mode.  Above  80  Km 
the  presence  of  the  shorter-wavelength  (2,4)  mode  is  evident, 
due  primarily  to  mode  coupling  processes  involving  the  mean 
zonal  wind  and  meridional  temperature  gradient  distributions. 

Over  most  of  the  0-140  Km  altitude  regime  the  largest 
density  perturbations  generally  occur  at  equatorial  latitudes 
for  the  diurnal  and  semidiurnal  tides.  However, 
considerable  latitude  structure  is  exhibited  at  any  given 
level.  A  discussion  was  also  included  of  the  possible  effect 
of  non-migrating  tidal  components,  (i.e.,  those  which  do  not 
migrate  with  the  apparent  motion  of  the  sun).  It  was  noted 
that  wind  data  from  the  Jicamarca,  Peru  (12°s)  incoherent 
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scatter  radar  indicate  the  presence  of  non-migrating 
components  with  amplitudes  as  much  as  an  order  of  magnitude 
larger  than  the  migrating  contribution.  Thus  the  potential 
exists  (at  least  at  some  longitudes)  of  significantly  larger 
tidal  density  variations  in  the  low-latitude  stratosphere 
than  given  for  the  migrating  components.  It  was  recommended 
that  SHAD  density  data  from  the  HEAO  experiment,  and  in  the 
possession  of  the  MIT  Draper  Laboratories,  should  be  able  to 
reveal  non-migrating  tidal  components  in  stratospheric 
density  if  they  exist. 

The  density  work  was  subsequently  expanded  upon  to 
provide  inputs  to  the  point-analysis  and  Air  Force  Reference 
Atmosphere  efforts,  and  was  written  up  in  as  a  proposed  AFGL 
report  entitled  "Estimates  of  Point-Density  Errors  due  to 
Atmospheric  Tides  between  70  and  120  Km"  by  J.M.  Forbes. 
This  latter  report  includes  various  contours  of  "percent 
density  departures  from  mean  values  due  to  tides"  as  a 
function  of  height,  latitude,  local  time  and  season. 

3.  TEMPERATURES,  PRESSURES,  AND  DENSITIES  IN  THE  80-120  Km 

REGION 

3.1  Introduction 

At  the  present  time  the  standard  reference  for 
temperature,  density,  and  wind  specification  between  80  and 
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120  Km  is  the  1972  COSPAR  International  Reference  Atmosphere 
(CIRA  1972).  Temperatures  and  densities  are  given  as 
monthly  averages  at  5  Km  height  increments  every  10°  in 
latitude  for  the  Northern  Hemisphere.  However,  data  are 
extremely  sparse  in  this  region;  as  Groves  (see  CIRA  1972) 
points  out,  CIRA  1972  suffers  from  the  same  limitations  as 
CIRA  1965  above  60  Km; 

(i)  Data  from  all  longitudes  are  combined  without 
consideration  of  longitudinal  effects.  Most  data 
are  from  N.  American  sites,  and  so  any  longitudinal 
bias  would  be  towards  the  W.  Hemisphere. 

(ii)  Insufficient  S.  Hemisphere  data  were  available  for 
developing  a  separate  model.  Therefore,  S. 
Hemisphere  data  were  combined  with  N.  Hemisphere 
data  with  a  six-month  change  of  date. 

(iii)  Due  to  insufficient  data,  no  account  was  taken  of 

local  time  in  development  of  the  models. 

Consequently  the  temperature  and  density  fields  may 
be  diurnally-biased  as  well. 

The  current  Air  Force  Reference  Atmosphere  (Cole  and  Kantor, 
1978)  only  extends  to  90Km.  Development  of  the  temperature 
and  density  inputs  between  80  and  120  Km  for  the  new  Air 
Force  Reference  Atmosphere  (AFRA)  has  involved  building  upon 
the  existing  CIRA  1972  model  data  base  (see  section  3.2). 
Since  the  greatest  abundance  of  experimental  data  is  in  the 
form  of  temperatures,  and  since  temperature  is  one  of  the 
meteorologial  fields  (besides  winds)  for  which  we  have  a 
firmer  theoretical  and  intuitive  base  for  understanding  its 
behavior  and  structure,  the  temperature  field  is  the  basis 
for  development  of  the  model.  Following  the  procedure 


followed  in  CIRA,  given  a  specification  of  the  pressure 
field  at  some  lower  boundary,  the  density  field  at  higher 
altitudes  were  derived  from  the  barometric  and  ideal  gas 
laws.  Available  experimental  density  determinations  are  used 
to  provide  consistency  checks  on  these  densities.  A 
preliminary  analysis  of  this  data  was  written  up  in  the  form 
of  a  report  "Temperature  Structure  of  the  80  to  120Km  Region" 
by  J.M.  Forbes,  as  input  to  the  COSPAR  working  Group  on  the 
new  Cospar  International  Reference  Atmosphere.  The  report 
has  been  published  in  the  MAP  Handbook  Series,  Volume  XVI. 

3.2.  Description  of  the  Temperature  Data  Base 

The  following  is  a  description  of  the  data  which  has 
formed  the  basis  for  the  80-120  Km  specification  of 
temperature,  density  and  pressure  recommended  for  the  new 
AFRA: 

(i)  The  CIRA  1972  temperature  model  between  80  and  120 
Km  was  based  heavily  on  data  from  the  NASA 
Meteorological  Sounding  Rocket  Program  (MSRP) 
collected  prior  to  1967.  The  primary  techniques 
were  pitot  tube  and  grenade  measurements.  The  MSPP 
was  phased  out  in  1973  .  Between  1967  and  1973  32 
pitot  tube  and  135  rocket  grenade  experiments  were 
conducted  which  generally  yielded  temperature  data 
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above  80Km.  The  soundings  were  made  at  Wallops  I. 
( 38°N ,  75°W) ,  Ft.  Churchill  (59°N,  94°W) , 

Pt.  Barrow  (71°N,  157°W),  Natal  <6°S, 

35°W) ,  Arecibo  (18°N,  67°W) ,  Arenosillo 

(37°N),  Eglin  {30°,  87°W),  and  Kourou  (5°N, 

53°W)  (Smith  et _ al,  1967,  1968,  1969,  1970, 

1971) . 

(ii)  Gaigerov  et  al(1984)  present  extensive  analyses  of 
temperatures  from  Soviet  rocket  measurements 
(mostly  grenade  method)  north  and  south  of  the 
equator  in  the  Eastern  Hemisphere.  Most  of  the 
Southern  Hemisphere  data  were  collected  after  1970, 
and  hence  were  not  included  in  the  1972  CIRA. 
Gaigerov  et  al  also  discuss  the  results  of  several 
intercomparisons  with  other  independent 
measurements  of  temperature,  and  report  that 
adjustments  have  been  made  for  any  biases  which 
might  have  existed  in  their  raw  data.  Monthly 
average  data  between  80  and  90-100Km  are  available 
from  Heiss  I.  (81°N,  58°E) ,  Volgograd  (48°N, 

44°E),  Molodezhnaya  (68°S,  45°E),  and  from 

research  vessel  soundings  in  the  Pacific  near  0°, 
20°S,  40°S,  and  50°S.  Keeping  in  mind  known 

deficiencies  in  the  Soviet  data,  this  information 
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now  enables  us  to  examine  possible  longitudinal  and 
latitudinal  asymmetries  in  the  structure  of  the 
mesopause  region.  S.  Hemisphere  data  may  in  fact 
introduce  an  undesired  bias  in  the  AFRA  if  mixed 


with  N.  Hemisphere  data,  and  so  this  aspect  of  the 
analysis  has  been  pursued  with  caution. 

Five  years  (1970-1975)  of  temperature  profiles  from 
incoherent  scatter  measurements  in  the  E-region 
(100-130Km)  over  Arecibo,  Puerto  Rico  (18°N, 
67°W)  and  Millstone  Hill,  Massachusetts  (42°N, 
71°W)  were  available  for  analysis.  These  data 
have  been  analyzed  in  parts,  mostly  with  regard  to 
the  semidiurnal  oscillation  by  Salah  (1974),  Salah 
et  al  (1975),  Salah  and  Wand  (1974),  and  Wand 
(1976,  1983).  This  investigator  has  pooled  all 
the  available  data,  separated  mean  and  semidiurnal 
tidal  components,  and  constructed  monthly 
averages.  A  total  of  over  1,500  profiles  (each) 
are  available  from  Arecibo  and  Millstone  Hill. 
These  data  are  considered  extremely  important  in 
terms  of  "matching"  the  rocket-based  temperature 
structure  of  the  80-100Km  region  with  the 


A  list  of  stations  from  vhich  rocket  and  radar  data  has 
originated  for  the  new  AFRA  is  included  in  Table  I.  A 
preliminary  analysis  was  first  performed  to  ascertain  whether 
evidence  exists  for  longitudinal  and  latitudinal  asymmetries 
in  the  temperature  structure  of  the  80Km  to  lOOKm  region, 
and  whether  sufficient  data  are  available  to  delineate 
these  dependences  in  a  reference  atmosphere. 


!  TABLE  I 

l 

Locations  of  rocket*  measurements  and  incoherent  scatter  radar** 
[  measurements  which  form  the  basis  data  for  the  new  AFRA  between  80Km 

1  and  120Km. 


"Western  Hemisphere"  "Eastern  Hemisphere" 


V  Thule 

( 76°N,  69°W) 

Heiss  I . 

( 81°N,  58°E) 

Pt.  Barrow 

( 71°N,  157°W) 

Volgograd 

Sardinia 

(48°N,  44°E ) 

Ft.  Churchill 

(59°N,  94°W) 

{ 40°N,  10°E) 

'  **  Millstone  Hill 

(42°N,  71°W) 

Guam 

(13°N,145°E) 

•  Wallops  I. 

( 38°N,  75°W) 

Kwajalein 

(9°N,168°E) 

£  White  Sands 

(32°N,  106°W) 

Thumba 

( 8°N ,  77°E) 

r  Eglin 

( 30°N,  87°W) 

Res.  Vessels 

(0°) 

Cape  Kennedy 

( 28°N,  80°W) 

( 22°N,  159°W) 

Res.  Vessels 

(20°) 

Barking  Sands 

Carnavon 

(25°S,114°E) 

**  Arecibo 

( 18°N,  67°W) 

Woomera 

( 31°S , 1 36°E ) 

1  Antigua 

( 17°N,  62°W) 

Res.  Vessels 

( 50°S ) 

•  Ft.  Sherman 

( 9°N,  80°W) 

Kerguelen  I . 

(49°S) 

!  Kourou 

( 5°N,  53°W) 

Res.  Vessels 

(50°) 

Natal 

Ascension  I. 

(6°S,  35°W) 

(8°S,  14°W) 

Molodezhnaya 

( 68°S ,  45°E) 

1  *  With  some  exceptions,  data 

are  generally 

available  between 

80-100Km 


**  Data  generally  available  between  100-130Km 
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In  Figure  3  variations  in  monthly  temperatures  at  85Km 
for  individual  stations  are  compared.  The  comparison  between 
Pt.  Barrow  (71°N,  157°w)  and  Heiss  I.  (81°N,  58°E) 

suggests  5-10°K  higher  temperatures  at  Heiss  I.  in  the 
summer  and  5-10°  cooler  temperatures  between  January  and 

March.  Although  these  two  stations  are  separated  by  10° 

in  latitude,  the  discrepancy  is  opposite  to  what  one  might 

expect  from  the  positive  (negative)  pole-to-equator 

temperature  gradient  assumed  to  exist  in  Northern  Hemisphere 
summer  (winter)  months.  An  examination  of  vertical 
structures  at  the  two  stations  indicates  that  the  summer 
mesopause  minimum  is  near  90Km  at  Heiss  I.  as  opposed  to  85Km 
at  Pt.  Barrow,  and  this  is  in  itself  is  an  important 
contributor  to  their  differences  in  monthly  behavior  at  85Km. 
Although  the  Volgograd  (48oN>  44oE)  data  during 

summer  exhibit  10-15°K  higher  temperatures  that  Ft. 

Churchill  (59°N,  94°W) ,  their  11°  separation  in 

latitude  is  sufficient  to  account  for  more  than  half  of  this 
difference  assuming  a  realistic  latitude  gradient  in 
temperature  (see  following  figures).  Figure  3  also  suggests 
a  much  larger  temperature  gradient  in  the  eastern  than 

western  hemisphere,  but  it  must  be  remembered  that  Ft. 

Churchill  and  Pt.  Barrow  are  only  separated  by  12° 
latitude,  whereas  the  separation  between  Heiss  I.  and 

Volgorad  is  33o  latitude.  Obviously,  to  make  any 


convincing  statements  about  latitude  structure  we  should 
examine  all  possible  data  at  a  given  height.  This  will  be 
done  below.  Before  leaving  Figure  3,  note  that  Southern 
Hemisphere  data  at  40°S  and  50°  agree  quite  well  with 
the  Northern  Hemisphere  data  at  similar  latitudes. 

In  Figures  4  and  5  the  latitude  structures  of 
temperature  during  summer  (mostly  July)  and  winter  (mostly 
January),  respectively,  are  depicted  at  80,  90,  and  lOOKm. 
An  obvious  feature  of  these  plots  is  that  the 
Eastern/Western/Northern/Southern  Hemisphere  data 
collectively  delineate  fairly  well-defined  patterns. 

Figure  6  depicts  a  selection  of  temperature  measurements 
between  100  and  130Km  at  various  latitudes  during  July  and 
January.  Comparisons  are  made  with  the  MSIS-83  model  (Hedin 
et  al,  1983),  as  this  is  a  likely  candidate  for  the  new  AFRA 
above  lOOKm  or  so.  The  data  shown  in  Figure  6  during  January 
are  extremely  consistent  with  each  other  and  with  the  MSIS-83 
model,  and  do  not  exhibit  any  significant  latitude 
structure.  During  July,  however,  there  appears  to  exist  a 
significant  positive  equator-to-pole  temperature  gradient. 
At  115Km,  the  temperature  varies  from  about  268°K  at 
Kwajalein  (9°N)  to  320°K  at  Arecibo  (18°N)  to  370°K 
at  Wallops  I.  (38°N)  and  Millstone  Hill  (42°N).  A  small 
temperature  difference  (20°K)  of  this  sense  between  18°N 
and  42°N  is  specified  in  the  MSIS-83  model.  Compared  with 
the  earlier  AFRA,  the  present  data  set  provides  improved 
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Figure  3.  Temperature  vs.  month  at  85Km  for  various  stations  which 
allow  examination  of  possible  longitudinal  or  hemispherical  asymmetries 
within  specific  latitude  belts.  CIRA  1972  values  are  shown  for 
comparison. 


seasonal,  latitudinal,  and  longitudinal  coverage  in  the 
80-100Km  region,  and  a  combination  of  incoherent  scatter  and 
rocket  data  in  the  100-120Km  region  which  allows  a  much 
improved  delineation  of  lower  thermosphere  temperature 
structure.  At  the  same  time,  although  some  individual 
station  comparisons  indicate  measurable  asymmetries  in 
longitude  and  latitude,  data  are  still  insufficient  to 
separate  these  effects;  that  is,  to  provide  a  reliable 
description  of  latitude  structure  as  a  function  of 
longitude,  or  of  longitude  structure  at  any  given  latitude. 
It  is  therefore  recommended  that  in  order  to  obtain  the  best 
description  of  Northern  Hemisphere  temperatures,  pressures, 
and  densities  between  80  and  120Km,  one  must  combine  together 
data  from  all  longitudes  and  from  the  Southern  Hemisphere 
(with  a  6-month  change  of  date)  in  order  to  obtain  even  a 
minimally  acceptable  distribution  of  data.  Thus  the  AFRA 
model  input  as  it  is  provided  here  is  not  capable  of 
distinguishing  longitudinal  or  hemispheric  dependences  if 
they  exist.  Soviet  data  (which  has  been  bias-corrected)  has 
also  been  included  in  the  data  set.  Inclusion  of  Southern 
Hemisphere  and/or  Soviet  data  may  arouse  some  concern.  As 
it  turns  out,  both  of  these  data  show  remarkable  consistency 
with  the  North  American  data,  and  in  some  cases  fill  crucial 
gaps  or  extend  the  data  set  and  model  credibility  to  higher 
latitudes  than  would  otherwise  be  unattainable.  High 
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Figure  4.  Temperature  vs.  latitude  for  measurements  representative  of 
July  north  of  the  equator  in  the  Western  (•)  and  Eastern  (  ■*  > 
Hemispheres,  and  January  south  of  the  equator  inthe  Eastern  Hemisphere 
(4*).  Where  data  under  these  exact  conditions  were  no  available  in  the 
Western  Hemisphere,  data  points  form  8°  August  (O),  6°  August  (&), 
August  (♦),  and  June  (O)  were  inserted  to  allow  a  more  complete 
delineation  of  the  latitude  structure.  CIRA  1972  values  are  shown  for 
comparison. 


Figure  5.  Temperature  vs.  latitude  for  measurements  representative  of 
January  north  of  the  equator  in  the  Western  (•)  and  Eastern  (*) 
Hemispheres.  Where  data  under  these  exact  conditions  were  not 
available  in  the  Western  Hemisphere,  data  points  from  6°S  February  (O 
),  6°  December  (A),  and  February  (O)  were  inserted  toa  llow  a  more 
complete  delineation  of  the  latitude  structure.  CIRA  1972  values  are 
shown  for  comparison. 
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Figure  6.  Vertical  structures  of  temperature  from  lOOKm  to  130Km  at 
various  latitudes.  MSIS-83  values  are  shown  for  comparison. 
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latitudes  show  the  greatest  degree  of  variability,  so  it  was 
important  that  every  possible  piece  of  information  be 
considered  in  development  of  the  model. 

3.3  Method  of  Analysis 

The  data  set  (consisting  of  monthly  averages  at  the 
observing  sites)  is  distributed  unevenly  in  latitude,  height, 

and  by  month.  In  order  to  construct  monthly  tables  at  15° 
latitude  increment,  it  is  necessary  to  fit  the  data  in  a 

least-squares  sense  using  appropriate  functions,  and  then  to 

use  these  formulas  to  specify  the  data  at  15o  latitude 

increments.  For  a  given  month  and  altitude  the  best 
latitudinal  fit  with  the  least  number  of  terms  was  found  to 
be  given  by  a  function  with  the  following  form  (  0  * 
colatitude) : 

T(z)  *  T0 ( z )  +  Aicos20  +  A2cos30  + 

A3cos40 

The  fit  was  made  to  the  data  set  from  70  to  120Km  every  5Km 
for  every  month  of  the  year.  (A  total  of  11  X  12  *  121 
fits).  The  polynomials  were  then  used  to  generate 
temperature  tables  every  15°  latitude,  every  5Km 
(70-120Km),  for  every  month.  These  data  do  not,  however, 
represent  a  smooth  profile  at  any  given  latitude.  The  next 
step  was  to  fit  an  analytic  function  to  every  15°  latitude 
data  set  (for  every  month)  between  70Km  or  so  and  120Km  such 


31 


.  *  4  '  fc 


that  the  temperature  and  temperature  derivative  match  the 
corresponding  values  from  the  existing  USAF  Reference 
Atmosphere  at  some  lower  boundary.  specifically,  the 
following  procedure  was  found  to  work  well:  Only  temperature 
data  from  80  to  120Km  in  5Km  increments  was  included  in  the 
fit.  Additional  data  points  at  68,  72,  74,  and  78Km  were 
obtained  from  the  old  AFRA  tables.  A  fourth-order  polynomial 
in  altitude  was  fit  to  all  these  data  points.  In  all  cases 
the  fit  was  extremely  close  to  the  data  points  in  the  70-90Km 
region,  and  formed  a  well-defined  minimum  (the  mesopause) 
near  90Km.  Tables  were  then  generated  every  2Km  from  the 
polynomial  fit.  A  final  adjustment  at  the  lower  boundary 
was  made  such  that  the  68Km  point  (generated  from  the 
polynomial)  was  replaced  by  the  actual  temperature  value  from 
the  old  AFRA  tables,  and  the  value  at  70Km  was  replaced  with 
the  value  obtained  by  linearly  interpolating  between  the 
value  at  68Km  and  the  polynomial-generated  value  at  72Km. 
This  ensured  exact  matching  of  temperature  with  the  old  AFRA 
at  68Km,  very  close  (but  approximate)  matching  with 
temperatures  and  gradients  from  the  old  AFRA  between  70  and 
80Km,  merging  into  the  newly-defined  mesopause  region 
(85-95Km)  and  lower-thermosphere  (100-120Km)  regions  above. 

It  is  important  to  note  two  points  relating  to  the 
method  analysis  as  adopted  above.  First,  no  attempt  was  made 
to  impose  or  fit  any  functional  dependence  vs.  month  at  a 


given  height  or  latitude.  Second,  since  the  temperature 
field  at  120Km  specified  by  the  upper  portion  of  the  AFRA 
(which  will  probably  be  the  MSIS-83  model)  is  dynamic  and 
dependent  on  solar,  geomagnetic,  and  other  variations,  a 
single  match  could  not  be  performed  at  these  levels  to 
conform  with  the  tabular  nature  of  the  model  below  120Km. 

3.4  Tables  of  Temperature,  Pressure,  and  Density 

Tables  of  monthly  Temperature,  Pressure,  and  Density 
extending  from  70Km  to  120Km  in  2Km  height  increments  at 
0°,  15°N,  30°N,  45°N,  60°N,  and  75°N  are 
appended  to  this  report.  These  tables  are  in  exactly  the  same 
format  as  the  existing  Air  Force  Reference  Atmosphere  Tables 
(Cole  and  Kantor,  1978).  The  pressure  tables  were  calculated 
using  the  barometric  law  and  a  reference  pressure  at  68Km 
equal  to  that  given  by  the  AFRA  (Cole  and  Kantor,  1978).  The 
latitude  dependence  of  the  acceleration  due  to  gravity  at  sea 
level  was  taken  from  Cole  and  Kantor  (1978),  and  the  formula 
for  the  height  dependence  of  g  was  taken  from  the  1966  U.S. 
Standard  Atmosphere  Supplements.  The  height  dependence  of 
the  mean  molecular  weight  above  90  Km  was  adopted  from  the 
MSIS  model  (Hedin  et  al ,  1983)  for  equinox  conditions,  an 
F10.7  cm  solar  flux  of  100  and  Ap=15;  the  values  for  M  from 
80  to  120Km  in  steps  of  lOKm  were  as  follows:  28.96,  28.90, 


28.46,  27.42,  26.23  in  a.m.u..  A  linear  interpolation  was 

utilized  at  intermediate  altitudes. 

3.5  Variability 

In  addition  to  the  marked  seasonal  variability  depicted 
in  the  temperature,  pressure,  and  density  tables,  within  a 
given  month  or  even  a  single  day  there  exist  significant 
variations  about  mean  values  due  to  gravity  waves,  tides,  and 
perhaps  other  sources.  Yet,  while  the  Tables  are  based  on 
temperature  because  it  is  the  more  globally  available 
measured  field,  the  variability  characteristics  of  total  mass 
density  are  of  much  greater  Air  Force  operational  concern. 
Therefore,  density  information  from  several  typical  stations 
are  used  here  to  quantitatively  characterize  the 
variability  of  density  between  70  and  llOKm.  Specifically, 
the  data  employed  originate  from  NASA  MSRP  data  from  Pt. 
Barrow  (71°N),  Ft.  Churchill  (59°N),  Wallops  I.  (38°N) 
and  the  Natal  (6°S)  -  Ascension  (8°s)  pair.  (Smith  et 

al,  1964;  1966,1967,1968,1969,1970,1971;  Theon  et  al,  1972). 
The  monthly  and  local  time  distributions  from  each  of  these 
sites  is  given  in  Figure  7.  Of  the  227  total  soundings 
represented,  207  are  grenade  soundings  (which  yield 
temperature  and  horizontal  winds  as  the  primary  variable, 
from  which  pressure  and  density  are  derived)  and  20  pitot 


probe  soundings  (which  actually  measure  density).  While 
there  exist  scattered  falling-sphere  density  data  which  will 
enable  important  consistency  checks  on  the  current  density 
tabulations,  the  NASA  MSRP  data  represent  the  best 
statistical  base  to  provide  information  on  density 
variability  as  a  function  of  season  and  latitude. 

Percent  Standard  deviations  of  temperature  and  density 
from  seasonal-mean  values  at  Pt.  Barrow,  Ft.  Churchill,  and 
Wallops  I.  are  given  in  Tables  II-IV.  Seasonal  means  are 
utilized  to  provide  better  statistics  given  the  small  number 
of  soundings  available.  Values  based  on  the  combined 
stations  of  Natal  (6°S)  and  Ascension  I  (8°S)  are  based 
on  annual  means  in  Table  V  due  to  the  small  seasonal 
variations  present  near  the  equator  (Theon  et  al ,  1972). 
Note  that  the  standard  deviations  for  temperature  are 
generally  smaller  than  for  density,  the  latter  values 
typically  lying  between  5-20%  between  70  and  llOKm. 

Theoretical  considerations  must  be  invoked  to  explain 
the  origin  of  the  above  variations.  Note  from  Figure  7  that 
in  toto  the  data  are  distributed  over  all  parts  of  the  day, 
with  some  bias  towards  the  1200-2000  hour  local  time  period. 
However,  within  a  given  season  the  local  time  distribution  is 
much  more  uneven,  essentially  precluding  determination  of 
tidal  components;  furthermore,  by  the  same  token,  the 
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tabulated  standard  deviations  do  not  fully  account  for  tidal 
effects.  Since  tides  represent  global  oscillations  with 
fairly  well-known  excitation  and  propagation  characteristics, 
useful  theoretical  estimates  of  this  variability  source  can 
be  obtained.  The  model  of  Forbes  ( 1982a, b)  has  been  utilized 
to  produce  such  estimates.  Figure  8  illustrates  percent 
density  departures  from  diurnal  mean  values  for  80  and  lOOKm 
between  0  and  60°  latitude.  An  important  feature  to  note 
is  the  small-amplitude  tidal  oscillation  predicated 
polevardof  about  54°  latitude.  This  implies  that  the 
density  variations  listed  in  the  TABLES  for  Pt.  Barrow 
(71°N)  and  Ft.  Churchill  (59°N)  may  be  presumed  to  be 
predominantly  a  reflection  of  gravity  wave  effects.  In  fact, 
the  height  dependence  of  density  standard  deviations  at  Pt. 
Barrow  is  similar  to  that  of  radar  echoes  (and  hence 
turbulent  irregularities)  observed  with  the  MST  Radar  at 
Poker  Flat  (65°N; 


nu  zi 

xstaarr  standard  dotations  of  seasonal  -  kan 

TEMEBRA3URES  AND  DENSITIES  AT  PT.  BARRON  (71%) 


ALT 

(Km) 

TEMP 
st.  a 

DENSITY 

ST.D. 

NO. 

GBSV. 

TEMP 

ST.D. 

SMSSBB 

DENSITY 

ST.D. 

NO. 

(BSV. 

TEMP 

ST.D. 

mnaam 

DENSITY 

ST.D. 

i 

NO. 

(BSV). 

70 

6 

20 

18 

3 

6 

10 

9 

24 

16 

74 

7 

19 

18 

2 

7 

9 

5 

26 

16 

78 

9 

17 

18 

2 

8 

9 

6 

29 

16 

82 

10 

13 

17 

2 

7 

8 

9 

29 

16 

86 

9 

12 

17 

4 

8 

8 

13 

26 

13 

90 

9 

10 

15 

7 

10 

7 

14 

21 

11 

94 

5 

7 

6 

24 

14 

2 

2 

11 

3 

37 


TALE  III 


I 


mean  stmtowd  db/iaticns  of  seasonal  -  tew 

TEMKRATJRES  AND  DENSITIES  AT  FT.  (BURCHILL  (59^) 


ALT 

(Km) 

TEMP 

ST.D. 

JQ22EEB 

DENSITY 

ST.D. 

NO. 

(BSV. 

TEMP 

ST.D. 

SMBB 

DENSITY 

ST.D. 

NO. 

(BSV. 

TEMP 

ST.D. 

SHttMG/FALL 
DENSITY  NO. 
ST.D.  (BSV. 

70 

5 

17 

29 

2 

6 

12 

3 

19 

13 

74 

7 

17 

29 

3 

6 

12 

3 

20 

13 

78 

8 

16 

28 

4 

7 

12 

5 

21 

13 

82 

9 

14 

27 

6 

9 

12 

6 

19 

13 

86 

6 

14 

26 

6 

9 

12 

6 

8 

12 

90 

11 

13 

19 

10 

15 

7 

12 

15 

9 

94 

5 

10 

11 

.5 

8 

2 

10 

3 

4 

98 

3 

9 

4 

- 

- 

1 

102 

7 

13 

4 

- 

- 

1 

106 

6 

7 

4 

- 

- 

1 

110 

16 

16 

4 

- 

- 

1 

114 

16 

8 

4 

- 

— 

1 

U8 

15 

11 

4 

- 

- 

1 

ABLE  3V 

mean  stanem©  dh/iatioe  of  seasonal  -  ban 

TEMERA3DRE  MID  DENSITIES  AT  WALLOPS  I.  (38<ta) 


WTWTER 

SOBER 

spruc/fail 

ALT 

TEMP 

DENSITY 

NO. 

TEMP 

DENSITY 

NO. 

TEMP 

DENSITY 

NO. 

(Km) 

ST.D. 

ST.D. 

(BSV. 

ST.D. 

st.  a 

(BSV. 

ST.D. 

ST.D.  (BSV, 

70 

4 

9 

28 

4 

9 

23 

4 

13 

38 

74 

7 

10 

25 

5 

8 

23 

6 

13 

38 

78 

7 

8 

21 

7 

9 

20 

7 

14 

36 

82 

7 

11 

16 

6 

9 

18 

8 

12 

32 

86 

6 

10 

10 

8 

14 

16 

8 

10 

31 

90 

11 

18 

6 

8 

13 

13 

14 

14 

19 

94 

- 

- 

1 

6 

8 

5 

3 

6 

3 

98 

9 

1 

4 

3 

.6 

2 

102 

14 

5 

3 

7 

1 

2 

106 

8 

19 

3 

2 

6 

2 

no 

11 

27 

3 

.1 

10 

2 

114 

6 

25 

2 

4 

11 

2 

U8 

5 

27 

2 

8 

11 

2 

1 


TfBLB  V 

PERCENT  STANDARD  D0/IAT  ICNS  OF  NCUAL  -  WAN 
TEM3RA3DRES  JWD  DENSITIES  AT  NATAL  -  A60NSIDN  (6  TO  8%) 


AT 
< Km ) 

TEMP 
ST.  D. 

DENSITY 
ST.  D. 

NO. 

CBS?. 

70 

4 

6 

32 

74 

5 

9 

31 

78 

6 

7 

30 

82 

5 

10 

30 

86 

8 

10 

30 

90 

11 

12 

29 

94 

14 

13 

12 

98 

22 

16 

3 

102 

23 

.1 

2 

106 

- 

- 

1 

Balsey  et  al,  1983);  that  is,  summer  is  characterized  by  the 
predominant  occurence  of  echoes  between  80  and  lOOKm,  where 
as  during 

winter  the  echoes  appear  to  originate  predominantly  below 
80Km.  The  turbulence  arises  from  convective  instabilities 
associated  with  gravity  waves  whose  level  of  breakdown  is 
modulated  by  the  seasonal ly-dependent  mean  winds  of  the 
troposphere  and  stratosphere. 

The  latitude  dependence  of  gravity-wave  related 
perturbations  in  the  70-100Km  height  region  is  not  known. 
However,  there  is  no  reason  to  expect  that  gravity  wave 
activity  should  increase  significantly  with  decreasing 
latitude.  The  data  in  Figure  8,  on  the  other  hand,  predict 
that  tidal  effects  increase  markedly  as  one  gets  closer  to 
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the  equator,  suggesting  that  tides  predominate  equatorward  of 
some  intermediate  latitude.  The  standard  deviations  of 
temperature  and  density  in  the  Tables  are  in  fact  consistent 
with  the  temperature  computations  presented  by 

Forbest 1982a, b)  and  the  densities  in  Figure  8.  At  middle 
latitudes,  say  30-55°  latitude,  it  is  not  unreasonable  to 
assume  that  tides  and  gravity  waves  contribute  about  equally 
to  the  observed  variability  of  the  region. 

In  the  100-120Km  regime  gravity  waves  most  likely  exist 
with  amplitudes  no  greater  than  their  convectively-limited 
values  at  lower  levels,  whereas  tides  experience  exponential 
growth  with  height  throughout  the  region.  In  addition, 
other  sources  of  tides  such  as  UV  and  EUV  solar  radiation 

absorption  assume  importance  above  95-100Km.  Estimated 
point  density  errors  at  120Km  as  presented  in  Figure  9 
commonly  occur  in  the  15-30  %  range.  This  order  of  magnitude 
is  supported  by  the  semidiurnal  temperature  perturbations 
obtained  from  an  analysis  of  Millstone  Hill  (42°N)  and 

Arecibo  (18°N)  temperatures  performed  by  Dr.  Forbes  and 
presented  in  Figure  10.  Figure  10  illustrates  seasonal 
averages  of  the  semidiurnal  temperature  oscillations  as 
determined  from  25  days  of  simultaneous  measurements  at 
Millstone  Hill  and  Arecibo.  These  data  represent  a  subset  of 
the  profiles  utilized  in  constructing  the  monthly  and 

diurnal-mean  profiles  discussed  in  Section  3.2.  The 


Millstone  Hill  amplitude  profile  is  characterized  by  a  peak 
near  110-115Km  which  varies  in  amplitude  from  45K  (summer, 
equinox)  to  60K  (winter).  The  Arecibo  temperatures 
generally  do  not  exhibit  a  well-defined  peak.  Both  stations 
exhibit  a  downward  phase  progression  with  characteristic 
vertical  scales  of  30-45Km.  Further,  the  Millstone  Hill 
phases  lead  those  at  Arecibo  by  about  4  hours  in  summer,  and 
lag  the  Arecibo  phases  by  4-6  hours  during  winter.  These 
amplitude  and  phase  characteristics  suggest  that  the 
observations  reflect  the  strong  presence  of  (2,4)  and  (2,5) 
semidurnal  tidal  modes  propagating  upwards  from  below  lOOKm. 

3.6  Discussion  of  Densities 

Since  total  mass  densities  represent  derived  quantities 
in  the  appended  tables,  it  is  necessary  to  perform 
consistency  checks  for  validation  purposes,  and  to  ascertain 
what  new  information  these  results  provide.  The  average 
density  profiles  corresponding  to  the  standard  deviations 
listed  in  Tables  III  and  IV  are  used  to  accomplish  this.  In 
addition,  since  MSIS-83  is  a  likely  candidate  for  the  AFRA  at 
satellite  altitudes,  comparisons  with  this  model  are  also 
made  here  to  aid  in  the  formulation  of  a  means  to  merge  the 
dynamic  model  above  and  the  tables  at  lower  heights  (see  also 
Section  3.4).  These  intercomparisons  are  performed  at  90, 
110,  and  120Km  in  Figures  11,  12,  and  13,  respectively. 
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Figure  10  Seasonal  averages  of  semidiurnal 


temperatures 


derived  from  25  simultaneous  daytime  experiments  at 
Arecibo  and  Millstone  Hill. 


The  present  model  is  seen  to  agree  extremely  well  with 
the  AFRA  at  90Km  (Figure  11).  Both  models  indicate  a 
predominantly  annual  variation,  with  a  greater  amplitude  at 
60°N  than  30°N.  The  only  marked  difference  is  that  at 
30°N  the  AFRA  predicts  the  density  maximum  to  occur  between 
February  and  April,  whereas  the  present  analysis  indicates 
maximum  densities  during  the  April-June  period;  both  models 
indicate  a  minimum  in  density  to  occur  during  the 
August-October  period.  At  60°  latitude  the  density  maxima 
for  both  models  occur  during  April-May  with  minima  during 
December-January .  The  annual  variation  predicted  by  MSIS-83 
is  in  good  agreement  at  60°N,  but  at  30°N  MSIS-83 
contains  little  variation  with  month  of  the  year.  The 
densities  at  Pt.  Barrow  appear  to  contain  a  significant 
semiannual  component,  and  about  half  the  variability 
indicated  by  the  other  sources  at  60°N.  The  Pt.  Barrow 
data,  on  the  other  hand,  appear  to  confirm  the  annual 
variation  with  amplitude  and  mean  value  about  20%  less  than 
the  other  models  examined. 

At  HOKm  (Figure  12),  the  current  model  indicates  a  more 
or  less  annual  variation  at  high  latitudes,  but  with  maximum 
during  December-Apri 1  and  minimum  during  June-July  ...  an 
approximately  180°  phase  difference  from  the  annual  cycle 
at  90Km!  The  data  at  Ft.  Churchill  are  in  agreement  with  the 
phase  and  amplitude  of  this  variation,  but  with  a  slightly 
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Figure  11.  Data  comparisons  at  90  Km. 
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Figure  12.  Data  compari 


Pigure  13.  Data  comparisons  at  120  Km. 
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smaller  mean  value.  The  MSIS-83  model  also  agrees  with  the 
sense  of  this  variation,  except  with  a  significantly  smaller 
amplitude.  The  density  variation  at  30°  latitude  now 
reflects  a  significant  semiannual  component.  This  result  is 
not  inconsistent  with  data  at  Wallops  I.,  except  that  during 
winter  months  there  does  not  exist  sufficient  data  at  this 
altitude  to  confirm  the  expected  decrease  in  density.  The 
MSIS-83  variation  is  still  annual  at  30°N,  with  only  a 
10-12%  amplitude  and  30%  smaller  mean  value. 

At  120Km  the  semiannual  variation  at  30°  latitude  is 
particularly  predominant,  while  at  60°  latitude  the 
distortion  of  the  annual  variation  also  indicates  presence  of 
a  significant  semiannual  component  as  well.  The  Ft.  Churchill 
data  and  MSIS-83  model  yield  similar  overall  magnitudes 
similar  to  the  current  results  at  60°N,  but  are 
approximately  30%  less  at  30°  latitude. 

3.7  Merging  with  MSIS-83  above  120  Km 

The  above  results  suggest  that  MSIS-83  is  in  reasonably 
good  agreement  with  available  data  at  90Km,  but  differs 
measurably  at  HOKm  and  120Km,  particularly  at  30°  latitude 
where  20-30%  underestimates  are  evident.  Practical  use  of  the 
new  AFRA  requires  a  functionally  efficient  algorithm  for 


merging  the  lower-altitude  tables  (<120Km)  with  the  MSIS-83 
model  above.  In  addition,  the  MSIS-83  contains  other 
geophysical  variations  such  as  those  connected  with  solar  and 
geomagnetic  activity  which  are  not  relevant  at  lower 
altitudes.  It  is  recommended  that,  at  least  initially,  the 
merging  of  these  two  model  segments  be  accomplished  as 
follows: 

(1)  Store  tabular  data  below  HOKm 

(2)  Use  MSIS-83  computer  algorithm  down  to  130Km 

(3)  Between  HOKm  and  130Km  determine  density  and 
pressure  by  linearly  interpolating  the  natural 
logarithms  of  these  fields,  and  determine 
temperature  by  direct  linear  interpolation. 

The  above  method  will  naturally  damp  out  solar  -  geomagnetic 
variations  from  130Km  to  HOKm.  However,  internal  consistency 
between  pressure,  density,  and  temperature  will  not 
necessarily  be  maintained  between  HOKm  and  130Km.  This 
should  not  be  cause  for  concern  for  operational  utilization 
of  the  model.  In  addition,  the  first  derivatives  of  these 
fields  will  not  match  at  HOKm  and  130Km.  After  initial 
experimentation  with  the  merging  method  described  above,  more 
sophisticated  methods  of  smoothing,  or  modification  of  the 
exact  altitudes  of  matching,  may  warrant  examination. 


Unless  specific  solar,  geomagnetic,  and  geophysical 
inputs  to  MSIS-83  are  given,  it  is  not  possible  to  provide  a 
usable  or  useful  set  of  tables  covering  the  80-200km  regime, 
unless  these  would  be  constructed  solely  for  illustrative 
purposes;  this  is  because  the  matching  conditions  would  have 
to  be  changed  for  each  set  of  solar  geophysical  conditions 
specified.  Rather,  the  model  must  by  given  in  computer 
algorithm  form  for  those  who  would  ultimately  make 
operational  use  of  the  model.  For  these  reasons,  the  tables 
given  here  only  extend  to  120Km,  and  a  recommended  method  of 
merging  has  been  provided  for  implementation  in  the  working 
version  of  the  new  AFRA. 
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